Bogoliubov Excitations in a Kronig-Penney Potential 
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With use of the Kronig-Penney model, we study the excitation spectrum of a Bose-Einstein 
condensate in a one- dimensional periodic potential. We solve the Bogoliubov equations analytically 
and obtain the band structure of the excitation spectrum for arbitrary values of the lattice depth. 
We find that the excitation spectrum is gapless and linear at low energies, and that it is due to the 
anomalous tunneling of low energy excitations, predicted by Kagan et at. 
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I. INTRODUCTION 

In studies of Bose-Einstein condensation in ultra- 
cold atomic gases, elementary excitations play an im- 
portant role in understanding various properties of the 
Bose-Einstein condensates, such as dynamics, thermo- 
dynamics and superfluidity. Elementary excitations 
have been experimentally investigated in trapped Bose- 
Einstein condensates, including the observation of col- 
lective modes 0, 0, and the measurement of the Bo- 
goliubov excitation spectrum with use of Bragg spec- 
troscopy [J El- 

An optical lattice is a periodic potential for atoms cre- 
ated by standing waves of laser beams. Experimental 
observations have also revealed properties of elementary 
excitations of Bose-Einstein condensates in an optical lat- 
tice. Stoferle et al. have studied excitation spectra of 
Bose-Einstein condensates in an optical lattice 5j. The 
Bose-Einstein condensates were excited by a modulation 
of the lattice depth, and a broad continuum of the ex- 
citation spectrum was observed; they suggested that the 
broad continuum characterizes the superfluid phase. 

In the present paper, we study elementary excitations 
of Bose-Einstein condensates in an optical lattice with 
use of a Kronig-Penney potential. In previous theoreti- 
cal papers ISSlllillSIIllllllllIl, elementary ex- 
citations in a sinusoidal lattice potential have been stud- 
ied. Berg-S0rensen et al. numerically solved the Bogoli- 
ubov equations for a one-dimensional optical lattice, and 
they also calculated the low energy excitations analyti- 
cally for shallow lattices within the Thomas- Fermi limit; 
they showed that the excitation spectrum is phonon-likc 
at low energies Q . The phonon dispersion of the excita- 
tion spectrum is directly connected to the superfluidity of 
a Bose-Einstein condensate in an optical lattice |l5j , and 
most of the papers support the phonon dispersion. In 
contrast, Ichioka et al. have reported on softening of the 



'Electronic address: 
t Electronic address: 



danshita@kh.phys.wascda.ac.jp 



3tuchiyaQiis.u-tokyo.ac.jpr 



low energy excitations for deep lattices by numerically 
solving the Bogoliubov equations [l^. Such softening of 
the excitation spectrum, if it really exists, suggests an 
instability of superfluidity of a Bose-Einstein condensate 
in an optical lattice. One of our purposes is to investigate 
the possibility of the softening by analytically solving the 
Bogoliubov equations with a Kronig-Penney potential. 
We will show that the softening does not occur in the 
Kronig-Penney model. 

In previous papers, analytical methods to calculate the 
excitation spectrum in the presence of a sinusoidal lat- 
tice potential are limited cither to deep or shallow lat- 
tices. For deep lattices where the overlap of the con- 
densate wave functions in neighboring sites is sufficiently 
small, a tight-binding approximation provides analyti- 
cal expressions for the first band of the excitation spec- 
trum [1. 1^. Holim . For shallow lattices, the perturbative 
treatment of the lattice potential is applicable [l2L Il3| . 
An advantage of the Kronig-Penney model is that one 
can calculate the excitation spectrum exactly with ana- 
lytic methods for arbitrary values of the lattice depth. 

Another advantage of the Kronig-Penney model is that 
the excitation spectrum can be related to tunneling prop- 
erties of the excitations through a single potential barrier. 
Kagan et al. studied the tunneling problem of excitations 
and predicted that a barrier is transparent for excitations 
within a limited range of low energies; they called such 
a behavior the anomalous tunneling [lfj. We shall see, 
in fact, that the anomalous tunneling is crucial to the 
phonon-like form of the low energy excitations and that 
the softening of the excitation spectrum does not occur. 

The outline of the present paper is as follows. In Sec. 
II, we introduce a formulation of the problem using the 
Bogoliubov theory and calculate the condensate wave 
function in a Kronig-Penney potential. In Sec. Ill, we an- 
alytically solve the Bogoliubov equations and obtain the 
band structure of the excitation spectrum in the Kronig- 
Penney potential. We derive the phonon-like form of the 
low energy excitation. We also calculate the phonon ve- 
locity and the band gap. Conclusions and open questions 
are discussed in the final section. 
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II. CONDENSATE WAVE FUNCTION IN A 
KRONIG-PENNEY POTENTIAL 

We consider a Bose-Einstein condensate confined in a 
combined potential of an axisymmetric harmonic poten- 
tial and a one-dimensional periodic potential along the 
axial direction (the x axis). As the periodic potential, 
we adopt a Kronig-Penney potential, 



oo 

V{x) =V ^ S[x-na), 



(1) 



where a is the lattice constant, and Vq is the strength 
of the (5-function potential barrier, expressing the lattice 
depth. It is assumed that the condensate is so elongated 
along the axial direction that the axial confinement can 
be neglected. We assume that the frequency of the har- 
monic potential of the radial direction is large enough 
compared to the excitation energy for the axial direc- 
tion. In this situation, the one-dimensional treatment of 
the problem is justified. 

Our formulation of the problem is based on the mean- 
field theory which consists of the time- independent 
Gross-Pitaevskii equation and the Bogoliubov equations. 
They are 
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[i + V(x) + 2g\y {x)\ 2 . (4) 



-H 
2m dx 2 

Here fi is the chemical potential and m is the mass of 
an atom. Since the radial confinement is harmonic and 
sufficiently tight, the coupling constant is affected by the 
harmonic oscillator length a± of the radial confinement as 
g = 2fi %" , where a s is the s-wave scattering length [l7|. 
The Gross-Pitaevskii equation determines the conden- 
sate wave function ^>o(x). The Bogoliubov equations 
determine the energy e of the elementary excitations of 
the condensate and the wave functions of the excitation 
(u(x), v(x)) . The elementary excitations correspond to 
small fluctuations of the condensate wave function from 
the static equilibrium |lj ■ 

The periodic potential is sinusoidal in experiments of 
atomic gases trapped in an optical lattice. However, the 
Kronig-Penney model is useful to understand the prob- 
lem of the periodic potential qualitatively, because it al- 
lows an analytical treatment of the problem. Schematic 
picture of the condensate in the Kronig-Penney potential 
is shown in Fig.^ 

The chemical potential is related to the number of con- 
densate atoms No in each well by the normalization con- 
dition: 

"(n+l)a 



(5) 



At first, we shall solve Eq. J5J and obtain the con- 
densate wave function in a Kronig-Penney potential. We 
assume that the condensate does not have supercurrent; 
therefore one can regard, without loss of generality, the 
condensate wave function as real. Multiplying Eq. J5J by 



dx 



one obtains the first integral of Eq. J5J , 



dx 
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where £ = ^— is the healing length, and A is a con- 
stant corresponding to the value of the condensate wave 
function at the center of each lattice site. It is expressed 



A=^o{{n + \)a). 
Integrating Eq. (jSJ again, one obtains 



(7) 
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The boundary conditions at x = na and x = (n + 
determine the constants A and Xq- This solution has 
been obtained in Refs. [TslIT^]. 

In order to calculate the excitation spectrum analyti- 
cally, we assume that the lattice constant is sufficiently 
larger than the healing length. This assumption also al- 
lows us to relate the band structure of the excitation 
spectrum to the tunneling properties of the excitations. 
In this situation, the condensate wave function near the 
center of each lattice site is not affected by other po- 
tential barriers. Then, one approximately obtains the 
ground state solution of Eq. as 



T 1 \ rP-j. i f\x-na\+ x 
®a{x) = W^tanh M ^ 

n-^j a < x < ( n +^ ) " 



(9) 



where xq is determined by the boundary condition at 
x = na 



^o(na + 0) = ^ {na - 0), 



(10) 
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FIG. 1: The schematic picture of a condensate in a Kronig 
Penney potential. 
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FIG. 2: The solid line represents the chemical potential fi 
for a — 10£o, as a function of Vo. The dashed (dotted) line 
represents an approximate value of /_t when Vo S> fffto£o (Vo -C 
gn £ ). 
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This condensate wave function corresponds to the A 
limit of the Eq. ©. 

Substituting Eq. © into the normalization condition 
of Eq. (J5J, we derive the relation between the chemical 
potential and Nq: 



i- 2 i 



2— tanh— - 



gno, 



(13) 



where no = ^ is the averaged density of the condensate. 
One can obtain approximate solutions of Eq. I|13|l in the 
limits of Vb > gn £, and V < gn £ , where £ = ^g^ - 
When Vb 3> <?no£o, we expand Eq. 1|13|1 into power series 
of & an d ffiofe d bt in 

a Vo ' 



v-gn 
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In a similar way, when Vb <C gno^o, we expand Eq. i|13|) 
into power series of — and _Y° f , and obtain 
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In Eqs. (|14|l and l)15|l. we express the expansions up to 
the third order of the small parameters. We show the 
chemical potential as a function of the potential strength 
Vb in Fig. [2J The chemical potential increases mono- 
tonically as the potential strength increases, because the 
presence of the potential barriers makes the effect of re- 
pulsive interaction more significant. 



III. EXCITATION SPECTRUM 

In this section, we shall solve the Bogoliubov equa- 
tions with the condensate wave function of Eq. @ and 
calculate the excitation spectrum of the condensate in a 
Kronig-Penney potential. 



A. Single barrier problem 

A periodic potential can be viewed as a periodic ar- 
ray of potential barriers. The band structure of a sin- 
gle particle in a periodic potential can be expressed in 
terms of the tunneling properties of the particle in the 
presence of a single barrier potential [2(j. When one 
solves the Schrodinger equation with a single barrier, the 
wave function including the transmission amplitude is 
obtained. By imposing the Bloch's theorem on the wave 
function, one obtains an equation to determine the band 
structure [2fj|. 

We shall determine the band structure of the excitation 
spectrum. In a similar way to the case of a single particle, 
one can relate the band structure of the excitation spec- 
trum to the tunneling properties of the excitations in the 
presence of a single potential barrier, as discussed in the 
next subsection. In fact, once one solves the Bogoliubov 
equations with a single barrier and obtains the tunneling 
properties of the excitations, the band structure of the 
excitation spectrum can be easily obtained. Hence, we 
first focus our attention on the |a;| < | region, and solve 
the Bogoliubov equations analytically, with regard to the 
single barrier problem. 

Kagan et al. have solved the Bogoliubov equations with 
a single rectangular barrier and numerically calculated 
the transmission amplitude 0] . Although they have also 
approximately obtained an analytical expression for the 
transmission amplitude, the approximation is not valid 
at very low energies. Since behaviors of the transmission 
amplitude around zero energy are crucial to the form of 
the excitation spectrum in an optical lattice, we need to 
obtain another analytical expression for the transmission 
amplitude. In the model of a 5-function potential barrier 
which corresponds to the thin barrier limit of a rectangu- 
lar barrier, one can calculate the transmission amplitude 
exactly. 

There exist two independent solutions of the single bar- 
rier problem, corresponding to the types of scattering 
process. One solution ip l (x) describes the process where 
a Bogoliubov excitation comes from left, and the other 
solution ip r (x) describes the process where a Bogoliubov 
excitation comes from right. The schematic pictures of 
the solutions of the single barrier problem are shown in 
Fig. [3J Substituting the condensate wave function of Eq. 
101 into the Bogoliubov equations, let us obtain the solu- 
tions ijj l (x) and tp r (x). On both sides of the barrier, one 
obtains four particular solutions analytically [16| under 
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FIG. 3: The schematic pictures of the solutions of the single 
barrier problem. Scattering processes in which (a) the exci- 
tation comes from the left hand side, and (b) the excitation 
comes from the right hand side. Solid curves, arrows and the 
dashed curves represent the condensate wave functions, the 
scattering components and the localized components, respec- 
tively. 



the boundary condition at |x| 3> £: 

(ti„(x), Unfa))* °C 

These solutions are 



(16) 



u n (x) = A n e R 



. |x| + x \ . p n £, 
' 1 j ~ lSg ^ X >2he 



x I e + /i — fj, tanh 



2 f\ x\ +x 
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v n (x) = A„e « 



, . |x| + x \ . , 
I 1 I -isgn(x) — 



x I e — fj, + fj, tanh 
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where 



pi,2 = ±p=±V2m(vV+?-M), (19) 
P3,4 = T«7 = =RV 2m(vV + £ 2 + a»), (20) 



/2fi+sgn(x) i\/Ef 



, n = 1, 



A„ = ^ V^-^n(x)iy/E; n = 2 , (21) 



1, n = 3,4 



-Ed 



2m' 



(22) 



Wave functions (iti(x), i>i(x)) and (u 2 (x), i>2 (a;)) 
describe scattering components. Wave functions 



(1*3(2;), U3 (a;)) at x < and (114(1), ^(a;))* at a; > de- 
scribe the localized components around the potential bar- 
rier, because they decay exponentially at |x| 3> £• Wave 
functions (1*3(35), 1*3 (x))* at a; > and (114(2;), 1)4(2;))* at 
x < diverge far from the potential barrier. It is noted 
that Eqs. {T5)| and (|2*U)) can be derived by solving 



^(f^ + 2M). 
2m 2m 



(23) 



Equation (|23|) expresses the Bogoliubov spectrum for a 
uniform system. The normalization constant A„ of the 
scattering components is determined to satisfy 



u n {x), v n (x) 



E p + fi±e 



2s 



(24) 



at |x| > f. 

We omit the unphysical divergent components. The 
solutions ip l (x) and ip r ( x ) are expressed as superposi- 
tions of the remaining components. The solution rp (x) 
is written as 



^(x) 



where the coefficients r, b, t, and c are the amplitudes 
of the reflected, the left localized, the transmitted, and 
the right localized components, respectively. They are 
functions of the energy e and the potential strength V . 
The boundary conditions at x — yield four equations 
to determine all the coefficients 




lf>'(+0) = V'(-0), 



dip 1 
dx 



+0 



dip 1 
dx 



2m V 1 



h 2 



no). 



(26) 
(27) 



These equations are linear simultaneous equations for the 
coefficients r, 6, t and c, and one can analytically solve 
them. Since the exact solutions of Eqs. (|26|1 and l|27|) are 
unnecessarily complicated, we only write approximate 
forms of the coefficients. When e -C \i and Vo /li£, 
one approximately obtains all the coefficients as 



Z : 



b = 



t 

Z = eVo - + 



Z 



(28) 

(29) 

(30) 
(31) 



Thus, we obtained ip l (x) analytically, and we can also 
calculate ip r (x) in the same way. 
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the phase shift S approaches zero as the energy is re- 
duced to zero. This means that the potential barrier is 
transparent for low energy excitations. This behavior of 
the low energy excitations has been called the anomalous 
tunneling by Kagan et al. jl6j. Equations (|30l) and l|31() 
shows that the peak of T has a Lorentzian shape with half 
width Ae ~ Vq, and the anomalous tunneling becomes 
restricted to only the excitations with very low energies 
as strength of the potential barrier increases. This pe- 
culiar tunneling behavior of the excitations appears for 
arbitrary values of the potential strength. 



FIG. 4: The transmission coefficient \t\ with Vo = 3/i£ (dot- 
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FIG. 5: The phase shift S of t with V = 3(4 (dotted lines), 
Vo = 10/i£ (solid lines) and Vo = 50p£ (dashed lines), as 
functions of the energy. 



The transmission coefficient T = \t\ 2 and the phase 
shift S = arg(i) are shown in Figs. 01 and respectively, 
as a function of the energy. Expanding t around e = 0, 
one can analytically obtain approximate expressions of 
\t\ and 6 



1 - a 



(i 



The coefficients a and (3 are 



(32) 
(33) 



2(^o~K)(^o 3 + ^o 2 + 2K^) 2 ~4( M e) 3 ) + 9(^o) 2 
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where 



-(34) 



(35) 



(36) 



It is obvious from Figs. 0] and [5] and Eqs. (|32[l and l|33|) 

that the transmission coefficient T approaches unity and 



B. Band structure of excitation spectrum 

In this subsection, wc analytically calculate the band 
structure of the excitation spectrum with use of the solu- 
tion of the single barrier problem obtained in the previous 
subsection. We show that the tunneling properties of the 
excitations determine the band structure. 

Since the Bogoliubov equations arc linear differential 
equations, a general solution of the equations can be de- 
scribed as a linear combination of independent solutions 
with the same energy. We can write a general solution 
of the Bogoliubov equations in the region |x| < | as a 
linear combination of ip l (x) and ijj r (x) pi) : 



x < 



(37) 



where x an d C are arbitrary constants. 

Now we extend this solution to all regions of x by 
means of the Bloch's theorem. The Bloch's theorem as- 
serts that ip satisfies 



tp(x + a) = e h tp( x ): 



dx 



e R 



dx 



(38) 
(39) 



where q is the quasi-momentum. Equations l|38|l and l|39|) 

at x = — t| yields an equation expressing the relation be- 
tween the excitation energy e and the quasi-momentum 



(qa 

cos l¥ 



-T Z ipa 1 

e n + —e 

2t 2t 



(40) 



Since we are assuming o>(, the localized components 
around the potential barriers, which decay exponentially 
and vanish at \x\ — % , do not appear explicitly in Eq. 

The Wronskian defined as 

WW*,ij i )=u j *—u i -u i ^u j *+v j *^v i -v i -^v j *Ul) 
dx dx dx dx 

yields relations between r and t. We can simplify Eq. 
(|40|l by the relations. One can easily prove from the Bo- 
goliubov equations that W is independent of x when ip J 
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FIG. 6: The band structure of the excitation spectrum in 
the presence of a Kronig-Penney potential with a — 10£o and 
Vo = 5<?no£o are shown. Shaded areas express band gaps due 
to the absence of solutions of Eq. 11441 . (a) The band structure 
is shown, (b) The left-hand side of Eq. 14411 as a function of 
energy is shown. 



and ip l have the same energy. By evaluating W(if) , 
one obtains the conservation law of the energy flux [l 



\t\ 2 + \r\ 2 



1. 



(42) 



Meanwhile, by evaluating W(ip r *> w)i one obtains the 
relation 



t = \t\e i5 ,r = ±i\r\e iS . 



(43) 



Substituting Eqs. JUJ and |@3J into Eq. one obtains 
the simplified relation between the excitation energy and 
the quasi-momentum 



-r-. = cos — 

\t\ \ h 



(44) 



This relation is ex actly the same form as that in the case 
of a single particle [20| . We clearly see from Eq. (|44l) that 
the tunneling properties in the single barrier problem de- 
termine the relation between the excitation energy and 
the quasi-momentum, namely the band structure. 

Solving Eq. I|44(l for the excitation energy e, we obtain 
the band structure of the excitation spectrum as shown 
in Fig. Ela). Fig. HJb) shows the left-hand side of Eq. 
(|44|l as a function of e. There exists no solution of Eq. 
(|44|l when the absolute value of the left-hand side ex- 
ceeds unity, because the absolute value of the right-hand 
side is equal to or less than unity. Therefore, the energy 
regions where the absolute value of the left-hand side is 
greater than unity correspond to the forbidden regions of 
the excitation energy, namely the band gaps. They are 
expressed as the shaded regions in Fig.|SJ 

In the strong potential limit Vo — > oo, the widths of all 
the energy bands become narrow. The energy bands ap- 
proach the discrete eigenenergies of the Bogoliubov equa- 
tions for a single potential well, because the condensate 
is perfectly divided into each well. 
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FIG. 7: The phonon velocities of Eq. I|46^ in Kronig-Penney 
potentials with a = 10£o and 20£o, as a function of the po- 
tential strength Vo- 



C. Phonon dispersion 

Fig. Eta) shows that the first band of the excitation 
spectrum is phonon-likc in the low energy regions. We 
shall verify that the form of the excitation spectrum at 
low energies is phonon-like for arbitrary values of the po- 
tential depth. The excitation spectrum in the low energy 
limit e — > is calculated by expanding Eq. (|44|l around 
e = 0. Substituting p ~ Eqs. and into 

Eq. H44fl . one obtains the phonon dispersion of the exci- 
tation spectrum 



e ~ cq, 

where the phonon velocity c is 



(45) 



' Y m((a + /?£) 2 -2£ 2 a)' 
When Vq » <?no£o, Eq. iffify can be approximated as 



(46) 



Co 



gn Q a 



2V + gn a 



5^0 
2a 



(47) 



where c = [ s the phonon velocity in uniform sys- 

tems. When Vo <C gn ^o, Eq. (|4l)|> can be approximated 
as 



co 1 



3Vo 2 



16a£ (gn ) 



(48) 



The phonon velocities for a = 10£o and 20£o as func- 
tions of Vo are shown in Fig.0 The phonon velocity de- 
creases monotonically as the barrier strength increases, 
and this behavior is qualitatively consistent with the case 
of a sinusoidal potential [T(| . 

Thus the anomalous tunneling properties of Eqs. 
and H33|) . namely the perfect transmission of very low 
energy excitations, leads to the phonon dispersion. It 
is pointed out in Ref. |l4j that the excitation spectrum 
shows softening at large wavelengths for a deep sinusoidal 
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FIG. 8: The energy band gap between the lowest and the 
second lowest band for a — 10£o, as a function of the potential 
strength Vq. The dashed (dotted) line shows an approximate 
value of the band gap in the limit of Vb 3> gnoa (Vq -C gno^o)- 



lattice potential, relating to the instability of the super- 
fluidity. In contrast, the excitation spectrum is always 
phonon-like in the Kronig-Penney potential because the 
anomalous tunneling occurs even for very strong poten- 
tial barriers. The phonon dispersion of the excitation 
spectrum at low energies reflects the superfluidity of the 
condensate. In this sense, the anomalous tunneling is 
crucial to the stability of the superfluidity in the Kronig- 
Penney potential. 

However, our results cannot eliminate the possibility of 
the softening in a sinusoidal periodic potential, because 
we have assumed in our calculation that the lattice con- 
stant is sufficiently larger than the healing length. In the 
Kronig-Penney model, this assumption assures that the 
size of the condensate wave function in each well hardly 
changes as the lattice potential becomes deep; conse- 
quently the size of the condensate is always sufficiently 
larger than the healing length. In contrast, even when 
the lattice constant is sufficiently larger than the healing 
length, the size of the condensate wave function in each 
well can be comparable to or smaller than the healing 
length in a deep sinusoidal potential. This is because the 
deeper the sinusoidal lattice potential is, the tighter the 
confinement of each well is. Thus, a possibility of the 
softening is remaining in the situation where the size of 
the condensate in each well is comparable to or smaller 
than the healing length. 



D. Band gap 

The jth band gap Aj t h is determined by the excitation 
energies of the top of lower band and the bottom of higher 
band, and they are obtained by solving Eq. <|44|) at q = 
Here we show the first band gap Ai st between the 
lowest and the second lowest band, as a function of the 
potential strength Vq in Fig. [HI 

When Vo ^> gnga, the first band gap is much larger 
than the first band width. In this case, one can approxi- 



A 



1st. - 



2£c 



Vna 
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When Vq <C gno£o, one can obtain 



A, 



.(49) 



(50) 



Equation l|5U|l and Fig. [S] show that at first the band gap 
increases linearly with the potential strength. The first 
term of the right-hand side of Eq. 1(4*91 corresponds to the 
lowest excitation energy of a condensate in a single well. 
Since the second term vanishes in the limit of Vq — » oo, 
the band gap becomes the lowest excitation energy in a 
single well. 



E. Tight-binding limit 

When the potential is sufficiently strong Vq 3> gn Q a, 
one can approximately solve Eq. (|44|l and obtain the first 
band of the excitation spectrum. Since the excitation 
energy in the first band is much smaller than the chemical 
potential in this situation, the approximate expression of 
t of Eqs. I|3U|) and (|31|l is adequate. One can rewrite Eqs. 
(OP and (HI) as 



cos S = 1*1 f 1 



sin 8 = \t\- (—-2 



where 



t = 



^/e^V a (V -2^) + ^e 



(51) 
(52) 

(53) 



Substituting Eq. JHU, and (JHSJ into Eq. one 
obtains 



gn \ 



aVn 



1 



5^o 
2a 



sm 



(-) 

\2h) 



(54) 



where fi is replaced with gno by using Eq. (|14f) . 

Meanwhile, the first band of the excitation spectrum 
has been calculated using the tight-binding approxima- 
tion in many theoretical papers [H, 0, 10, 11J. For a very 
deep lattice, the spectrum takes the form 



sm 



(-) 

\2h) 



(55) 



where k is the compressibility of the system and 77 is the 
lowest energy bandwidth of stationary current-carrying 
states of the condensate ^(j- The compressibility k is 
defined by the thermodynamic relation 



1 

K 



'dNo~' 



(56) 
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The parameter r\ is related to the effective mass m* 
through the relation 



2Ti 2 



V 



The definition of the effective mass is 
f d 2 e I 



dk 2 



(57) 



(58) 



fe=0 



where k is the quasi-momentum of the current-carrying 
condensate and e is the energy per particle of the con- 
densate. We should be careful not to confuse the energy 
dispersion of the Bogoliubov excitations s(q), which is 
our major interest, and the energy per particle e(fc) of 
the condensate itself. 

Using Eqs. (Spl and (|5?jl . one can obtain 



f 

K 

V 



gn 



1 



a 



aVn 



1 



4?o 
a 



(59) 
(60) 



Here we have used Eq. I|f4ll to calculate the compress- 
ibility. Substituting Eqs. JHU and (JHOl into Eq. iJSSj), we 

see that our calculation of the first band of the excitation 
spectrum is consistent with the result of the tight-binding 
approximation when the potential is sufficiently strong. 



obtained, which determines the band structure. We have 
shown that the excitation spectrum is gapless and linear 
in the low energy region due to the anomalous tunneling 
property of the excitations. 

While we have considered current-free condensates 
in our calculations, the elementary excitations of con- 
densates with large current in periodic potentials are 
known to exhibit the Landau and dynamical instabili- 
ties 0, HE E3- Since both instabilities are associ- 
ated with the specific properties of the excitation spec- 
tra in the low energy region, the tunneling properties 
at low energy may affect the instability when the prob- 
lem is considered with the Kronig-Penney model. Hence, 
it will be interesting to study the excitations of the 
current-carrying condensates in the Kronig-Penney po- 
tential and discuss the relation between the instability 
and the anomalous tunneling of the Bogoliubov excita- 
tions. 

Note added. After the submission of this paper , w e 
learned of a relevant paper by Kovrizhin [24| • In Ref. |24| , 
the tunneling problem of the Bogoliubov excitations 
through a cfeZta-function potential barrier is investigated. 
Exact form of the wave function Eq. I)25|l is derived, and 
the anomalous behavior of the transmission coefficient is 
discussed. 
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